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Using first-principle simulations for the probability density of finding a '^He atom in the vicinity 
of the screw dislocation in solid ''He, we determine the binding energy to the dislocation nucleus 
Eb = 0.8 ± O.IK and the density of localized states at larger distances. The specific heat due to "^He 
features a peak similar to the one observed in recent experiments, and our model can also account 
for the observed increase in shear modulus at low temperature. We further discuss the role of "^He 
in the picture of superfluid defects. 
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PACS numbers: 67.80.-s, 05.30.Jp, 67.30.hm, 61.72.Lk 

The observation of a non-classical moment of ro- 
tational inertia (NCRI) in torsional oscillator experi- 
ments [1[ in solid *He has revived the debate on super- 
solids, but its understanding proved to be challenging. 
Several groups over the world have confirmed NCRI, but 
aspects such as the pressure dependence [ll H, disor- 
der [3|, y, history dependence (hysteresis) [5|, la], "crit- 
ical" velocity [l|, |5[, crystal growth [3|, oscillation fre- 
quency [5[, rim velocity, and ^He concentration [7| all 
show unexpected behavior and defy any simple physical 
picture (for a review, see [S, Isj). Theoretically, consen- 
sus is growing towards a network of superfluid defects 
as the mechanism of superflow, but some effects cannot 
yet be explained properly, especially at the quantitative 
level |]. 

One of the main puzzles is the effect of even minute ■^He 
concentrations. There is mounting evidence that the in- 
terplay between '^He impurities and crystallographic de- 
fects (dislocations in single crystals) is not innocuous and 
can, in fact, be understood theoretically to a large degree. 
In this Letter, we investigate this topic and focus on the 
■^He binding to screw dislocations, which are common 
crystal defects in solid ^He. 

Day and Beamish observed an increase in the shear 
modulus of the ^He crystal when the temperature is low- 
ered [10|, which could be understood from binding of '^He 
to dislocations at low temperatures. According to the 
Granato-Liicke theory [llj , dislocations move in response 
to shear stress in their glide plane. More precisely, they 
bow out between pinning centers provided by impuri- 
ties or intersections, which can reduce the shear modulus 
by 30% from its intrinsic value in a frequency indepen- 
dent way. When '^Hc binds to a dislocation, it acts as 
an additional pinning center. Since the change in shear 
modulus is quadratic in the length between the pinning 
centers, the shear modulus quickly recovers its intrinsic 
value. Remarkably, the shear modulus dependence on 
temperature is nearly identical to that of NCRI. Yet, the 
two phenomena are distinct: the NCRI signal can not be 
fully accounted for by the elasticity effect, nor can the 



reduction in NCRI by a factor of 100 when blocking an 
annulus be explained by elasticity arguments [l|, |9| . 

In the torsional oscillator experiments by Kim et al. 
[7| a minimum "^Hc concentration of the order of X3 ^ 
Ippb seems needed in order to observe NCRI. Then, 
NCRI increases until X3 '^ Ippm where a maximum is 
reached and finally disappears again for concentrations 
of about lOOppm. Specific heat measurements showed 
a nearly constant term in the specific heat at low tem- 
peratures scaling with the '^He concentration [I4I. After 
subtraction of the phonon contribution and the mysteri- 
ous constant term, a peak in the specific heat was found 
around T = 0.06K, which was claimed to be indepen- 
dent of X3 and indicative of the supersolid transition [l^l 
(see however [13|). In this Letter we show, however, that 
binding of '^He impurities to dislocations results in a spe- 
cific heat peak in the same temperature range. 

Our approach is numerical and based on Fcynman's 
path-integral formulation of quantum mechanics. The 
integrals over the paths (world lines) are efficiently eval- 
uated by the worm algorithm [14| . which has been suc- 
cessful in describing properties of crystallographic defects 
in solid *IIe [H, 0, [13 El • We now describe how the 
method needs to be modified in order to deal with '^He 
impurities. 
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expressed as a path integral with the usual discretization 
of the imaginary time (inverse temperature) (3 into M 
slices ((5 = /3/Af), 
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where R = {Rq, Ri, ■■., Rm = PRq) is a particular 
world-line configuration with Rj = {'■i.j,'"2.j, ...,^7V.j} 
the coordinates of all N particles in time slice j, and 
dR = dRo...dRM-i- AH permutations P of the bosons 
are incorporated in the periodic boundary condition 
Rm ~ PRo- We use the primitive approximation [19| 
where U{R) contains only the inter-particle interaction 
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FIG. 1: (Color online) Reweighting function (line) and mass 
histogram (bars) as a function of ratio between the impurity 
particle mass and the ''He mass. 



given by the Aziz potential [20[. The kinetic term T{R) 
is a product of free-particle propagators [l9| , 
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We cannot just add a substitutional "^He atom to the 
setup and wait for it to hop around, because the exchange 
amplitude between ^He and ^He atoms is very small in 
solid ^He, J34 ~ 10~^K [21| compared to the tempera- 
tures of interest (T ~ 0.5K). A partial solution to this 
problem is to allow for a special Monte Carlo update, 
which relabels '^He and ''He trajectories (a ^Hc trajec- 
tory which is not part of any exchange cycle is chosen 
at random) thus leaving the world line configuration un- 
changed. As both "^Hc and '^He interact via the same 
potential, only the kinetic part in Eq. ([T]) is affected by 
this update. Using Eq. ([2]), the acceptance probability is 
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= min{l, e-('^-'*)(™*-™3)/<5}^ (3) 



with li = J2i=o ('"fci.i-i-i ~ ^fei,i)^^^/2 and rrii the mass, 
where the index i = 3, 4 refers to the ^He and ^He particle 
of the current update, respectively. Typical acceptance 
ratios are of the order of 10~^ and thus prohibitively 
low. The problem is that the ^He trajectory has a bigger 
fluctuation volume than the ^He one. 

To overcome this problem we introduce an update 
which gradually changes the mass m of an impurity par- 
ticle over the interval [m^ , 7714] . We do not allow for more 
than one impurity atom and work with a discrete set of 
11 impurity masses m = m3 + Am*i where i = 0, 1, ... 10 
and Am = (7714 — m^)/10. A gradual change in mass al- 
lows the crystal to relax and readjust the configuration to 
the new impurity mass. If m is close or equal to mi, the 
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FIG. 2: (Color online) The probability density of finding a sin- 
gle He atom at a radial distance r (cylindrical coordinates) 
from the core of the screw dislocation for different tempera- 
tures. Errors are indicated by the dashed lines for T = 0.5A', 
and are of the same order at higher temperatures. The grey 
line shows the reference for ^He in the absence of an impurity. 
The density is 71 = 0.0295A~^. At large distances, r > 30A 
(r > 25 A for T > IK), the *He atoms are treated as inert 
particles (with fixed straight world lines) and form a zone 
inaccessible to the impurity atom. 



exchange updates are frequently accepted. The quanti- 
ties of interest are only measured in the "physical" sector, 
where m = 7713. The acceptance probability for changing 
the mass from 771 to tti ± Am, is 



Pm~.7n±A7n = min{l, (1 ± Am / m) 



3M/2 -l^Am/S 



}• (4) 



This in itself does not solve the problem, since on aver- 
age Pm->m-Am < 1 and we only rarely visit the low-mass 
sector of the configuration space. The final solution is in 
employing a reweighing (importance sampling) technique 
which ensures that the probability of visiting different 
mass sectors are approximately equal. This is achieved 
by introducing the reweighing function .9(771) shown in 
Fig. [1] into the acceptance probability 

Pm->m±Am — > Pm->rn±Arng{m) / g{m ± AtTi). (5) 

This enables the impurity to efficiently sweep over the 
entire mass range. Every time the impurity is a true 
•^He atom, we measure its distance r to the nucleus of 
the screw dislocations and update the histogram for the 
radial probability density g{r) shown in Fig. [2] 

Relating g{r) to the effective potential energy E{r) be- 
tween a "^He atom and the dislocation core is straightfor- 
ward. The exchange matrix element between ^He and 
''He is negligible (~ lO^^K) compared to the tempera- 
tures of interest (T > 20 mK) and can be ignored al- 
together, leaving us with the classical Boltzmann dis- 
tribution g{r) oc exp[— /3i?(r)]. At large distances we 
assume -^(7") ex r~^ from elasticity theory [22| and 
proceed as follows: first we we fit the tail of g{r) to 
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FIG. 3: (Color online) Specific heat as a function of temper- 
ature for a single ^He impurity in the strain field of a screw 
dislocation obtained directly from the g{r) curve at T = 0.5K 
in Fig.[2]and assuming Xd = lO^'/a^ = 7.6 x 10* cm"^. The 
specific heat is exponentially suppressed at temperatures well 
above where the maximum is reached. The dashed lines cor- 
respond to the errors coming from the error on g{r) in Fig. [5] 
The specific heat curves resulting from g{r) at higher temper- 
atures (T = IK and T = 1.5K in Fig. [2} are the same as the 
one shown, within error bars. 



goo exp(— /3i?/r^) law to determine the asymptotic be- 
havior (^oo and B are fit parameters), and then we ob- 
tain the potential energies from E{r) = —T\n[g{r)/goo] 
for distances r < 20A. The partition function is found 
by integrating exp[— /3i?(r)] over all lattice sites up to 
some cut-off value Tmax = l/^Trxd where Xd is the dislo- 
cation density per a^ and a the inter-particle distance. 
The specific heat cy shown in Fig. [Sj is directly cal- 
culated from g{r) shown in Fig. [2] with only Xd as a 
free parameter. The partition function can be written 
approximately as Z « l/xd + NBexp{PEB) where Nb 
is the number of the deepest binding sites (with energy 
—Eb ~ —Tln[gjnax/gcx.] = 0.8±0.1K) per lattice period. 

The specific heat maximum is roughly at Tmax ~ 
Eb/ l'n-[l/NBXd]) with only a logarithmic dependence on 
the free parameter Xd- The peak falls in the same tem- 
perature range as the peak of Ref. [12] . The peak ampli- 
tude scales with the '^He concentration (since the Pauli 
exclusion can be neglected for low x^ and assuming full 
equilibration) and has the typical shape of a Schottky 
peak for a system with two degenerate energy levels (see 
above). Further refinements of our model, such as work- 
ing with discrete lattice points close to the core, using a 
Fermi function, or a distribution of binding energies for 
different dislocation types, are possible, but do not seem 
needed. 

The Shevchenko state of the network of interconnected 

predicts a crossover in the 



superfluid dislocations [17|, \2c 
specific heat from cy ex T above the bulk transition tem- 
perature Tc to cy (X T^ at low temperatures T < Tc- 
This signal, however, might be extremely small and unde- 
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FIG. 4: The reduction in the winding numbers (W^ ) (pro- 
portional to the superfluid density) along the nucleus of the 
screw dislocation when the "^He concentration at the disloca- 
tion core is increased. The temperature is T = 0.5K and the 
density is n = 0.0295A~='. 



tectable leaving the '^He contribution as the leading one. 
Also, binding of '^He to dislocations has an immediate 
effect on the shear modulus leading to crystal stiffening 
at r < Tniax as observed in Ref. [1C| . 

Increasing the '^He concentration in the nucleus of the 
superfluid screw dislocation is expected to reduce the su- 
perfluid density. We assume that '^He atoms will cluster 
at the point where dislocations intersect. We model their 
effect by introducing different numbers of impurities to 
the dislocation core next to each other (in this simula- 
tion we do not employ any of the special updates and 
reweighing mentioned above and thus the "^He cluster al- 
ways remains in the core). We see in Fig. |4]that about 
four '^He atoms are required to suppress the superfluid re- 
sponse along the core. This mechanism may explain the 
reduction of the superfluid response for concentrations 
X3 > 1 ppm observed in Ref. Jl]. 

Our last consideration is about the kinetic relaxation 
of '^He atoms. So far we assumed thermodynamic equi- 
librium which is not necessarily the case. NMR mea- 
surements established that the tunneling motion of '^He 
atoms in *He crystals is characterized by the hopping 
amplitude J34 '^ 10""' K. Any strain field producing 
an energy level bias between the nearest neighbor sites 
^ w a[dE / dr) much larger than ZJ34, where z is the 
coordination number, will localize '^He atoms. To move 
around, impurities have to exchange energy with the en- 
vironment. At low temperatures, the leading mechanism 
is provided by the one-phonon coupling and leads to hop- 
Jl^(^T/Q%, where Qd is the Debye 



ping rates r~ 

temperature 2j| . It is clear from the value of the bind- 
ing energy that in the vicinity of the dislocation core the 
condition <^ ^ 2:734 is definitely satisfied. The slowest 
rate is for ^ ~ ZJ34, i.e. for ^He to cross the boundary 
between the band motion and localized states. One can 



see, that the corresponding relaxation time is of the or- 
der of years at low temperature leading to sample history 
dependent effects, as observed in experiments. 

In conclusion, we have studied numerically the binding 
of ^He to the screw dislocation from first principles. We 
find a binding energy of 0.8±0.1K in agreement with pub- 
lished estimates. The binding of '^He impurity atoms to 
dislocation cores at low temperature results in a specific 
heat peak in the same temperature interval as observed 
experimentally in Ref . 12, ll3| , and may radically change 
superfluid properties of the dislocation network even at 
minute ^He concentrations. Our data also provide quan- 
titative support to the mechanism proposed in Ref. [lO[ 
as an explanation for the crystal stiffening at low tem- 
perature. 
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